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Abstract: We survey the phenomenological constraints on abelian gauge bosons having 
masses in the MeV to multi-GeV mass range (using precision electroweak measurements, 
neutrino-electron and neutrino-nucleon scattering, electron and muon anomalous magnetic 
moments, upsilon decay, beam dump experiments, atomic parity violation, low-energy neutron 
scattering and primordial nucleosynthesis). We compute their implications for the three 
parameters that in general describe the low-energy properties of such bosons: their mass and 
their two possible types of dimensionless couplings (direct couplings to ordinary fermions and 
kinetic mixing with Standard Model hypercharge) . We argue that gauge bosons with very 
small couplings to ordinary fermions in this mass range are natural in string compactifications 
and are likely to be generic in theories for which the gravity scale is systematically smaller 
than the Planck mass ~ such as in extra-dimensional models - because of the necessity to 
suppress proton decay. Furthermore, because its couplings are weak, in the low-energy theory 
relevant to experiments at and below TeV scales the charge gauged by the new boson can 
appear to be broken, both by classical effects and by anomalies. In particular, if the new gauge 
charge appears to be anomalous, anomaly cancellation does not also require the introduction 
of new light fermions in the low-energy theory. Furthermore, the charge can appear to be 
conserved in the low-energy theory, despite the corresponding gauge boson having a mass. 
Our results reduce to those of other authors in the special cases where there is no kinetic 
mixing or there is no direct coupling to ordinary fermions, such as for recently proposed 
dark-matter scenarios. 
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1. Introduction and summary of results 

New particles need not have very large masses in order to have evaded discovery; they can 
also be quite light provided they couple weakly enough to the other particles we do see. This 
unremarkable observation has been reinforced by recent dark matter models, many of which 
introduce new particles at GeV or lower scales in order to provide dark-matter interpretations 
for various astrophysical anomalies [1]. This model-building exercise has emphasized how 
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comparatively small experimental efforts might close off a wide range of at-present allowed 
couplings and masses for putative new light particles [2, 3]. 

Light spin-one bosons 

Spin-one gauge bosons are particularly natural kinds of particles to seek at low energies, since 
(unlike most scalars) these can have light masses in a technically natural way. Furthermore, 
their couplings are reasonably restrictive, allowing only two kinds of dimensionless interactions 
with ordinary Standard Model particles: direct gauge couplings to ordinary matter and kinetic 
mixing [4] with Standard Model gauge bosons. Most extant surveys of constraints on particles 
of this type assume the existence of one or the other of these couplings, with older studies 
studying only direct gauge- fermion interactions [5, 6] and later studies (particularly for dark- 
matter motivated models) [7, 8, 9, 10, 11] usually allowing only kinetic mixing. 

In this paper we have both motivational and phenomenological goals. On the phenomeno- 
logical side, we analyze the constraints on new (abelian) gauge bosons, including both direct 
gauge-fermion couplings and gauge-boson kinetic mixing. In this way we include all of their 
dimensionless couplings, which (if all other things are equal) should dominate their behaviour 
at low energies. We can follow the interplay of these couplings with one another, and how 
this changes the bounds that can be inferred concerning the allowed parameter space. In 
particular we find in some cases (such as beam dump experiments) that bounds derived un- 
der the assumption of the absence of the other coupling can sometimes weaken, rather than 
strengthen, once the most general couplings are present. 

Our motivational goal in this paper is twofold. First, we argue that the existence of gauge 
bosons directly coupled to ordinary fermions is very likely to be a generic and robust property 
of any phenomenologically successful theory for which the gravity scale is much smaller than 
the GUT scale [12, 13, 14]. Next, we argue that these gauge bosons often very naturally have 
extremely weak gauge interactions within reasonable UV extensions of the low-energy theory, 
such as extra-dimensional models [15] and low-energy string vacua [16]. Besides motivating 
the otherwise potentially repulsive feature of having very small couplings, the smallness of 
these couplings (together with the low value for the fundamental gravity scale) also naturally 
tends to make the corresponding gauge bosons unusually light. 

The remainder of this paper is organized as follows. The rest of this section, §1, briefly 
summarizes the basic motivational arguments and phenomenological results. §2 then provides 
a more detailed theoretical background that motivates the sizes and kinds of couplings we 
consider, which may be skipped for those interested only in the bounds themselves. In §3 
we briefly summarize the basic properties of the new gauge boson, with details given in an 
Appendix. By diagonalizing all kinetic terms and masses we identify the physical combination 
of couplings that are bounded in the subsequent sections. The next three sections, §4, §5 and 
§6, then explore the bounds on these couplings that are most restrictive for successively lighter 
bosons, starting at the weak scale and working down to MeV scales. 
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Motivational summary 

Why consider light gauge bosons that couple directly to ordinary fermions? And why should 
their couplings be so small? We here briefly summarize the more lengthy motivations given 
below, in §2. 

Low-scale gravity and proton decay 

Weakly coupled gauge bosons are likely to be generic features of any (phenomenologically 
viable) UV physics for which the fundamental gravity scale is systematically small relative to 
the GUT scale, Mqut ~ 10^^ GeV. Such bosons arise because of the difficulty of reconciling 
a low gravity scale with the observed stability of the proton. After all, higher-dimension 
baryon- and lepton-violating interactions that generically cause proton decay are not ade- 
quately suppressed if they arise accompanied by a gravity scale that is much smaller than 
Mqut- Similarly, global symmetries cannot themselves stop proton decay if the present lore 
about the absence of global symmetries in quantum gravity [17] should prove to be true (as 
happens in string theory, in particular [18, 19]). 

This leaves low-energy gauge symmetries as the remaining generic mechanism for sup- 
pressing proton decay. Indeed, extra gauge bosons are often found in string vacua, and when 
the string scale is much smaller than the GUT scale, Ms <C MquTj these bosons typically play 
a crucial role in protecting protons from decaying. Furthermore, very weak gauge couplings 
appear naturally in such string compactifications, once modulus stabilization is included. In 
these systems the gauge couplings can be small because they are often inversely proportional 
to the volume of some higher-dimensional cycle, whose volume gets stabilized at very large 
values [20]. Similar things can also occur in non-stringy extra-dimensional models [21]. 

Unbroken gauge symmetry without unbroken gauge symmetry 

We believe there is a generic low-energy lesson to be drawn from how proton decay is avoided 
in phenomenological string constructions. This is because in these models, even though 
proton decay is forbidden by conservation of a gauged charge, the gauge boson that gauges 
this symmetry is not massless [16]. This combines the virtues of an unbroken symmetry 
(no proton decay), with the virtues of a broken symmetry (no new forces mediated by a 
massless gauge boson). ^ Usually this happy situation arises in the string examples because 
the gauge symmetry in question is anomalous, if judged solely by the light fermion content, 
with anomaly freedom restored through Green-Schwarz anomaly cancellation. But in four 
dimensions Green-Schwarz anomaly cancellation relies on the existence of a Goldstone boson, 
whose presence also ensures that the gauge boson acquires a nonzero mass. 

For these constructions the effective lagrangian obtained just below the string scale from 
matching to the stringy UV completion is invariant under the symmetry apart from an 
anomaly- cancelling term that breaks the symmetry in just the way required to cancel the 

'^Superconductors are similar in this regard: the photon acquires a mass without implying gross violations 
of charge conservation. 
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fermion loop anomalies. §2 argues that this property remains true (to all orders in pertur- 
bation theory) as one integrates out modes down to low energies. Leading symmetry break- 
ing contributions arise non-perturbatively, exponentially suppressed by the relevant gauge 
couplings. Consequently they remain negligibly small provided only that the gauge groups 
involved in the anomalies are weakly coupled. Although supersymmetry also plays a role in 
the explicit string examples usually examined, our point here is that this is not required for 
the basic mechanism that allows massive gauge bosons to coexist with conservation of the 
corresponding gauge charge. 

Phenomenological summary 

We next summarize, for convenience of reference, the combined bounds obtained from the 
constraints examined throughout the following sections. 

Mass vs coupling 

Fig. 1 presents a series of exclusion plots in the Ox — Mx plane, where ax = 5^/47r is the 
gauge-fermion coupling and Mx is the gauge boson mass. Each panel shows these bounds 
for different fixed values of the kinetic mixing parameter, sh?] (for details on the definition 
of variables, see §3). The figure shows the collective exclusion area of all of the different 
bounds considered in this paper. For concreteness they are calculated for a vector-like charge 
assignment, = X/a, with the choice X = B — L denoted by a lighter shading and the 
choice X = B denoted with a heavier shading. Comparison of the cases X = B and X = B—L 
shows how much the bounds strengthen once direct couplings to leptons are allowed. 

For = 0, the dominant bounds are from neutrino scattering, upsilon decay, anomalous 
magnetic moments, beam-dump experiments, neutron-nucleus scattering and nucleosynthesis. 
Once kinetic mixing is introduced, many of these bounds improve, with the exception of the 
beam-dump bounds. Once shr/ > 0.06, kinetic mixing becomes sufficiently strong that the 
VF-mass bound prevails over any other bounds in the Mx ^ Mz region. For shry = 1, we 
discard the region where the oblique T parameter is large (for details, see §4), and focus on 
the region where Mx > 385 GeV. In this region, it is the neutrino-electron scattering bound 
and the VF-mass bound that dominate. 

Mass vs Mixing angle 

It is useful to show these same bounds as exclusion plots in the mixing-angle/boson- mass 
plane, for fixed choices of the gauge-fermion coupling, ax- This allows contact to be made 
with similar bounds obtained in the context of dark matter-inspired U{1) models [9, 2, 3, 10], 
which correspond to the ax — )• limit of the bounds we find here. This version of the 
plots is shown in Figure 2, restricted to the MeV-GeV mass range (in order to facilitate the 
comparison with earlier work). 

For small, but non-zero, gauge coupling (ax ~ 10~^^) the bounds from beam dump 
experiments weaken significantly. However, another strong bound from neutrino-electron 
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M, (GeV) 

Figure 1: Summary of the constraints presented herein. Each plot shows the bound on the new gauge 
coupling, as a function of for various values of the kinetic-mixing parameter, shry, assuming a 
vector coupling Xf^^ — Xf^, ■= X, with X = B — L {X = B) drawn as sparse (dense) cross-hatching. 



scattering also begins to take effect. This bound dominates for larger ax, and once ax ^ W 
the entire MeV— GeV mass range is excluded. 
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Figure 2: Summary of the constraints on kinetic mixing relevant in the MeV-GeV mass range. Each 
plot shows the bound on the kinetic mixing parameter sh?] as a function of Mx, for Ux = 0, 1 x 10""'^'' 
and 1 X 10^^. The plot assumes a coupling Xi^ — X^^ = — 1, such as would be true \i X = B — L. 
Hatched regions are excluded. 

Since the bounds in Figure 2 all rely on coupling to leptons, in the case where X = B the 
constraints arise through the kinetic mixing and are independent of ax- The resulting plot 
for X = B \s therefore the same as is shown in the figure for ax = 0. However, as the gauge 
coupling is increased the neutron-nucleus scattering bound — discussed in §6.4 — eventually 
becomes important, first being visible as an exclusion in the sh.r} ~ Mx plane in the panel for 
ax ^ 10"^ in Fig. 2. 

2. Theoretical motivation 

This section elaborates the motivations for weakly coupled, very light gauge bosons alluded 
to above. This is done both by summarizing the consistency conditions they must satisfy 
within the low-energy effective theory relevant to experiments, and by describing how such 
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bosons actually arise from several representative UV completions in string theory and extra- 
dimensional models. 

2.1 Low^-energy gauge symmetries, consistency and anomaly cancellation 

Very general arguments [22, 23] indicate that the interplay between unitarity and Lorentz 
invariance require massless gauge bosons only to couple to conserved charges that generate 
exact symmetries of the matter action. Consequently we normally expect the direct couplings 
of very light gauge bosons to be similarly restricted. This section reviews these arguments, 
emphasizing how they can break down [24, 25, 26] if the energy scale. A, of any UV completion 
is sufficiently small compared with the gauge boson mass, M, and coupling, g: A < AiiM/g. 
For the present purposes it suffices to restrict our attention to abelian gauge bosons (see 
however [26] for some discussion of the nonabelian case). 

The upshot of the arguments summarized here is that massive spin-one bosons can couple 
in an essentially arbitrary way if their mass, M, lies within a factor g/A-K of the scale of UV 
completion. But once M becomes smaller than gK/AiT, then the corresponding boson must 
gauge an honest-to-God, linearly realized exact symmetry. In particular this symmetry must 
be anomaly free. However any anomalies that Standard Model fermions give a putative 
new gauge charge needn't be cancelled by adding new, exotic low-energy fermions; they can 
instead be cancelled by the Goldstone boson whose presence is in any case required if the 
gauge boson has a mass. But this latter sort of cancellation also requires the UV completion 
scale to satisfy A < AirM/g. 

Notice that for any given M the condition A < inAd/g need not require A to lie below 
the TeV scale if the coupling g is small enough. For instance, if M ~ 1 MeV then A lies 
above the TeV scale provided g < 10~^ (an upper limit often required in any case by the 
strong phenomenological bounds we find below). And, as subsequent sections argue, such 
small couplings can actually arise in a natural way from reasonable UV completions. 

Massless spin-one bosons 

What goes wrong if a spin-one particle is not coupled to matter by gauging an exact symmetry 
of the matter action? If the spin-one particle is massless, then the problem is that one 
must give up either Lorentz invariance or unitarity (provided the particle has non-derivative. 
Coulomb-like couplings that survive in the far infrared). Lorentz invariance and unitarity 
fight one another because the basic field, A^(x), cannot transform as a Lorentz 4- vector if Af^ 
creates and destroys massless spin-one particles [22, 23]. Instead it transforms as a 4- vector up 
to a gauge transformation, — )• A^ + d^uj, and so interactions must be kept gauge invariant 
in order to be Lorentz invariant [27]. 

Massive spin-one bosons 

For massive spin-one particles the argument proceeds differently, as is now described. The 
difference arises because a 4- vector field, A^, can represent a massive spin-one particle [22]. 
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To examine the relevance of symmetries, it is worth first considering couphng massive 
spin-one particles to other matter fields, ip, in some arbitrary non-gauge-invariant way, with 
lagrangian density C{Af^,ip). The first observation to make is that any such a lagrangian 
can be made gauge invariant for free, by introducing a Stiickelberg field, 0, according to the 
replacements Afj_ — t- Afj, := — d^cj) and ip ^ ^ := exp[—i(j)Q]'ip, where Q is a hermitian 
matrix acting on the fields ip. With this replacement the lagrangian C{Afj,, ^) is automatically 
invariant under the symmetry A^ — )■ -|- (J^w, (p ^ cp + co and ip — >• exp[za; Qjtp, since both 
Afj, and ^ are themselves invariant under these replacements. The original non-symmetric 
formulation corresponds to the specific gauge (p = 0. For gauge symmetry, absence of gauge 
invariance is evidently equivalent to nonlinearly realized gauge invariance (similar arguments 
can also be made in the nonabelian case [26]). 

But this gauge invariance is obtained at the expense of introducing a new scale. Since <p 
is dimensionless, its kinetic term involves a scale, v, 

Ain = ^M-^^' - y (d^^ - ^m) (d^^ - An . (2.1) 

In (p> = gauge the scale v is seen to be related to the gauge boson mass by the relation 
M = gv. In a general gauge the scale v controls the size of couplings between the canonically 
normalized field, ip = (pv, and other particles. For instance the coupling 

-Ccoupiing = -i(V^7^QV) - ^) , (2.2) 

shows that the {Tp^^Q'ip)d^'^ coupling is dimension-five, being suppressed by the scale v = 
M/g. 

Lagrangians with nonrenormalizable couplings like this must be interpreted as eff'ective 
field theories, whose predictive power relies on performing a low-energy expansion in powers 
of -E/A, for some UV scale A. The interpretation of the scale v then generically depends on 
the how high A is relative to AnM/g = Attv. We consider each case in turn. 

Light spin-one bosons: M <^ gK/Air 

If the gauge boson is very light compared with the UV scale, then its low-energy interactions 
should be describable by some renormalizable theory. But renormalizability is only consistent 
with a dimension-five interaction^ like the {ip'^^Qip)d^ip coupling of eq. (2.2) if this coupling 
is a redundant interaction, such as it would be if it could be removed by a field redefinition. 
A sufficient condition for an interaction of the form J^dfj_<p to be redundant in this way is if 
the field equations for ip were to imply the quantity J^{ip) satisfies d^J^ = Q [28]. This shows 
that if the gauge boson is to be arbitrarily light relative to A, its low-energy, renormalizable 
couplings must be to a (dimension-three) conserved current. This is the usual prescription 

^The careful reader will recognize that this argument assumes negligible anomalous dimensions, and so 
needs re-examination for strongly coupled theories. 
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for obtaining these couplings by gauging a linearly realized matter symmetry, for which J'^ is 
the usual Noether current. 

More generic massive spin-one bosons: M > gA/Air 

If, on the other hand, the dimension-five coupling {ip^^Qij))d^ip is not redundant, then there 
must be an upper bound on the UV scale: A <, 4ttv = AnM/g. Sometimes this may be seen 
from the energy-dependence predicted for the cross section of reactions in the low energy 
theory: if cr{E) oc l/(47r?;)^ then this would be larger than the unitarity bound a ^ l/E"^ 
for energies S ^ A ~ 47ru, indicating the failure at these energies of the low-energy approxi- 
mation. If so, the full UV completion must intervene at or below these energies to keep the 
theory unitary. 

The upshot is that spin-one particles can couple fairly arbitrarily to matter provided they 
are massive, and provided the energy scale. A, of any UV completion satisfies A < AnM/g, 
where M is the gauge boson mass and g is its coupling strength. (Everyday examples of 
spin-one particles of this type include the p meson or spin-one nuclei.) It is only spin-one 
particles with M < gA/in that must gauge linearly realized symmetries. 

2.2 Anomaly cancellation 

Any new gauge symmetry — henceforth denoted U{l)x — must be an exact symmetry 
(though possibly spontaneously broken), and in particular must be anomaly free. This is true 
regardless of whether the symmetry is the linearly realized symmetry of a light gauge boson, 
or the nonlinearly realized symmetry of a massive gauge boson. 

Of particular interest in this paper are models where the new symmetry acts on ordinary 
fermions, because a robust motivation for thinking about light gauge bosons is the avoidance 
of proton decay in models with a low gravity scale (more about which below). In this case 
these ordinary fermions usually contribute gauge anomalies for the new symmetry, and an 
important issue is how these anomalies are ultimately cancelled. The two main anomaly- 
cancellation scenarios then divide according to whether or not anomalies cancel among the 
SM fields themselves, or require the addition of new particles. 

Anomaly cancellation using only SM fields 

The simplest situation is where the new gauge symmetry is simply a linear combination of one 
or more of the SM's four classical global symmetries — baryon number B, electron number 
Lg, muon number and tau number Lr- In this situation there are only two independent 
combinations of these symmetries that are anomaly free^ [29], corresponding to arbitrary 
linear combinations of the anomaly- free symmetries Le — and L^ — Lr'. 

X = a{Le - L^) + h{L^ - Lr) . (2.3) 

Of course, evidence for neutrino oscillations [30] make it unlikely that these symmetries are 

unbroken in whatever replaces the Standard Model in our ultimate understanding of Nature. 

^Notice that B — L carries a Standard Model anomaly in the absence of sterile right-handed neutrinos (see 
below) . 
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Anomaly cancellation using the Green-Schwarz mechanism 

If more general combinations B, L^, and are to be gauged, it is necessary to introduce 
new particles that can cancel their Standard Model anomalies. For a new U{l)x symmetry 
the minimal way to do this is to add only the Goldstone boson, which must in any case be 
present if the corresponding gauge boson has a mass (as it typically must to avoid mediating 
a macroscopic, long-range new force, whose presence is strongly disfavoured by observations 
[31]). For a U{l)x symmetry this can always be done using the 4D version [32] of the Green- 
Schwarz mechanism [33]. Besides its intrinsic interest, this is a way of cancelling anomalies 
that actually arises from plausible UV physics, such as low energy string models. 

In principle, there are four types of new anomalies that can arise in 4D once the SM 
is supplemented by a new gauge symmetry, ^7(1)^^. These are proportional to Tr[XXX], 
Tr[XXy], Tr[Xyy] and Ty[XG"'G'^], where the trace is over all left-handed fermions and X 
denotes the new symmetry generator, Y is Standard Model hypercharge, and G'^ represents 
the generators of the Standard Model nonabelian gauge groups, SU{2)i^ x SU{3)c, as well 
as the generators of Lorentz transformations. In four dimensions CPT invariance implies 
the absence of pure gravitational anomalies, and anomaly cancellation within the Standard 
Model ensures the absence of anomalies of the form Tr[yyy] and Tt[G°'G^G'^]. 

It is always possible to redefine the new symmetry generator, V := X + (^Y, to remove 
one of the two mixed anomalies. For instance, Tr[yyy] = Tr[XXy] + 2(^Tr[Xyy] can be 
made to vanish by choosing C appropriately (provided Tr[Xyy] does not vanish). It suffices 
then to consider only the case of nonzero anomalies of the form Tr[yyy] and Tr:[VG°'G°'], 
where G"" now includes also the generator Y. The anomaly then can be written in the G"- 
and Lorentz-invariant form^ 



where T is the 'quantum action' (generator of IPI correlations), the symmetry parameter is 
normalized by (5X^ = S^o; and the coefficients, Cx, Ca and c^, are calculable. Here Fy = dV = 
Fx + C Fy is the gauge-boson field strength for the generator X -\- C,Y, while Fa is the same 
for the Standard Model gauge bosons and R is the gravitational curvature 2-form. 

Given the coefficients Cx, Ca and c^, here is how 4D Green-Schwarz anomaly cancellation 
works [32]. Consider the gauge kinetic lagrangian, including the Stiickelberg field (p, 



''a similar formulation can be made using the anomaly in its 'consistent' form, rather than the 'covariant' 
form used in the text. 




(2.4) 



■mv 



- i^^A^r - ^Tr[G^,Gr] - y (5^.^ - X>^){d,^ - X,) 

+4>{cx[Fx +CFr)A (Fx + CFy) + CaTr[G'^ A G"] - c^Tr[i? A i?]} . (2.5) 
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Here i2inv denotes those parts of the lagrangian that are invariant under all of the gauge 
symmetries that are not written explicitly. The second line is not invariant under gauge 
transformations because 4> is not; its variation precisely cancels the fermion anomaly, eq. (2.4). 

An important observation is that the anomaly cancelling term is dimension-five, and so 
is not renormalizable. For instance, in terms of the canonically normalized field, (p = <j)v, the 
first anomaly cancelling term is >Canom = {^/ f)Fx A + • • • , where / = v/cx- As before, 
this implies the existence of a UV-completion scale. A, above which the low-energy effective 
description breaks down [25]. For weakly coupled theories typically A < Attv ~ AirM/g ~ 
Ancxf marks the scale where the fields arise that are required to extend the Goldstone boson 
to a linear representation of the symmetry. 

Perhaps the most interesting feature of cancelling anomalies with the Green-Schwarz 
mechanism in this way is that the lagrangian remains invariant under the U{\) symmetry, 
apart from the anomaly-cancelling term. This is interesting because it means that the cor- 
responding charge still appears to be conserved in the low-energy theory, despite the gauge 
field being massive. This opens up interesting phenomenological possibilities for the gauging 
of symmetries like U{1)b and U{1)b-l-: which appear to be conserved in Nature but which 
are also ruled out as sources of the new long-range force that a massless gauge boson would 
imply. 

One might worry that arbitrary symmetry-breaking interactions might be generated by 
embedding the anomaly cancelling interactions (or the fermion triangle anomaly graph) into 
a quantum fluctuation. For instance \{ X = B, so the new gauge boson couples to baryon 
number, then why can't some complicated loop generate a Ai? = ±1 interaction, 0±i, that 
can mediate proton decay? After all, this can be U{1)b invariant if it arises multiplied by a 
factor e^**^, which carries baryon number A.B = =Fl. 

The difficulty with generating this kind of interaction is that it must involve (f) undiffer- 
entiated. But if we restrict £anom to constant (j) configurations, it becomes a total derivative. 
For constant 0, the dependence of observables on (j) is similar to the dependence of observables 
on the vacuum angle, 9. Consequently it arises at best only non-perturbatively, proportional 
at weak coupling to a power of ~ exp[— Stt^/^^], where g is the anomalous gauge coupling. 
As a result the only potentially dangerous contribution of this type comes from the mixed 
X-QCD-QCD anomaly, which can generate nontrivial (/)-dependence once we integrate down 
to scales < Aq^d- This is not dangerous in particular for the classical symmetries, Lg, 
and Lr, since these do not have mixed QCD anomalies [29]. 

Anomaly cancellation using new fermions 

More complicated possibilities for new gauge bosons emerge if new, light exotic fermions 
are allowed that also carry the new X charge (and so can also take part in the anomaly 
cancellation). We briefly describe some features involving such new exotic particles, although 
they do not play any role in our later phenomenological studies. 

The simplest example along these lines \s X = B — L, which is anomaly- free provided 
only that the SM spectrum is supplemented by three right-handed neutrinos (one for each 
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generation). Furthermore, conservation of L is consistent with all evidence for neutrino 
oscillations, although it would be ruled out should neutrinoless double-beta decay ever be 
witnessed. 

A practical way in which such new fermions can arise at TeV scales is if the UV theory 
at these scales is supersymmetric. In this case the plethora of new superpartners can change 
anomaly cancellation in one of two ways (or both). They can either directly contribute to 
the anomalies themselves, and possibly help anomalies cancel without recourse to the Green- 
Schwarz mechanism. Alternatively, they can modify the details of how the Green-Schwarz 
mechanism operates if the UV scale, u, associated with it is larger than the supersymmetry 
breaking scale, Ms^sy. 

In particular, supersymmetry typically relates the kinetic term for the Stiickelberg field, 
(2.5), with a Fayet-Iliopoulos term in the scalar potential [35], 



where r is a dynamical field whose vev acts as the low-energy Fayet-Iliopoulos parameter; 
the qi are the charges of the fields (jji under the U{\) in question. In string examples the field 
r corresponds to a modulus of the compactification, which controls the size of a cycle in the 
internal geometry on which some branes wrap. We note that the vanishing of the D-term is 
consistent with vanishing vevs of the charged fields if r = 0, i.e. the symmetry survives as 
an exact global symmetry when the cycle size vanishes (the singular locus). Small values of 
the vev are obtained if the cycle size is small. 

2.3 Motivations from UV physics 

The above summary outlines some of the theoretical constraints on coupling ordinary fermions 
to very light gauge bosons. This section shows how very small couplings can naturally appear 
in well- motivated ultraviolet physics, such as extra-dimensional models or string vacua. In 
particular, they often arise due to considerations of proton stability in constructions for which 
the gravity scale is small compared with the Planck scale, as we now explain. 

Proton decay in low-scale gravity models 

One of the surprises of the late 20th century was the discovery that the scale, Mg, of quantum 
gravity could be much smaller than the Planck scale, Mp = (SvrG)"^/^ ~ 10^^ GeV [12]. From 
the point of view of particle physics this possibility is remarkable for several reasons. Most 
obvious is the potential it allows for experimental detection if it should happen that Mg is in 
the vicinity of the TeV scale [13, 14]. 

But there is a potentially more wide-reaching consequence that Mg <^ Mp has for the low- 
energy sector: the suppression by powers of Mg/Mp it allows for otherwise UV-sensitive ra- 
diative corrections [36]. This suppression arises because the contribution of short-wavelength 
degrees of freedom can saturate at M^, allowing their effects to be suppressed by powers of 
the gravitational coupling. 




(2.6) 
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The most precise examples of this are provided by string theory, in the regime where the 
string scale is low, Mg := Mg <C Mp [12]. String theory makes the suppression of UV-sensitive 
contributions precise by providing an explicit stringy ultraviolet completion within which the 
effects of the full UV sector can be explored. Large-volume (LV) models [20] are particularly 
useful laboratories for these purposes, since these systematically exploit the expansion of 
inverse powers of the extra-dimensional volume (in string units), V := (Vol)/^^ ^ 1, and it 
is ultimately these kinds of powers that enforce the suppressions of interest since M^jMp oc 

Proton decay — that is, its experimental absence — turns out to impose a very general 
constraint on any fundamental theory of this type, with Mg <^ Mp. It does so because having 
Mg very small removes two of the standard ways of keeping the proton stable in specific 
models. On one hand quantum gravity, and string theory in particular [18, 19], seems to 
preclude the existence of global symmetries, and this forbids ensuring proton stability by 
simply using a conserved global charge (such as baryon number). 

If Mg is too small then it is also unlikely that such a symmetry simply emerges by accident 
for the lowest-dimension interactions in the low-energy effective theory. The problem in this 
case is that we know that generic higher-dimensional interactions, 

= E ^ ' (2-7) 

^ Mg' 

eventually do arise in the low-energy effective theory, such as the standard baryon-number 
violating 4-quark operators arising at dimension = 6 [37] in the low-energy limit of grand- 
unified theories (GUTs) [38, 39]. But a dimension-six interaction of the form O/M^ generically 
contributes a proton-decay rate of order F ~ m^/M^, where rup is the proton mass, which is 
too large to agree with observations once M falls below Mquj. ~ 10^^ GeV. 

The way theories with Mg <C Mp usually evade proton decay is through the appearance of 
a gauged C/(l), whose conservation forbids the decay. Of course, to be useful the gauged U (1) 
that appears must couple to the proton or its decay products in order to forbid its decay. But 
because this means ordinary particles couple to the new gauge boson, it potentially introduces 
other phenomenological issues. If the gauge symmetry is conserved, why isn't the gauge boson 
massless? If the gauge boson is light, why isn't the new boson seen in low-energy observations? 
If the gauge boson is heavy, the corresponding symmetry must be badly broken and so how 
can it help with proton decay? Interestingly, extant models can naturally address both of 
these issues, and often the low-energy mechanism that is used is Green-Schwarz anomaly 
cancellation with gauge boson mass generated through the Stiickelberg mechanism described 
above. Sometimes this mechanism is also combined with supersymmetry to suppress the 
dangerous decays. 

The existence of these gauge bosons, their properties, and the way they evade the above 
issues, may be among the few generic low-energy consequences of viable theories with a low 
gravity scale: Mg <^ Mqut- 

Sample symmetries: 
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The simplest proposals for new low-energy gauge groups that forbid proton decay are either 
baryon or lepton number, X = B or X = L. If the anomalies for these symmetries due to 
Standard Model fermions are cancelled through the Green-Schwarz mechanism, then no new 
light particles are required besides the massive gauge boson itself. 

More complicated examples are possible if the low-energy theory at TeV scales is super- 
symmetric. In this case symmetries like B — L, that in themselves cannot forbid proton decay, 
can help suppress proton decay if taken together with supersymmetry [16]. (For instance, 
the parity R = (— that is usually used to suppress proton decay in the MSSM is a 
combination of fermion number and B — L.) 

More general combinations of B and L can also suppress proton decay in supersymmetric 
theories. Ref. [16] provides a list of the kinds of symmetries of this type that can be relevant 
to proton decay, as well as the conditions they must satisfy in order to have their anomalies 
be cancelled through the Green-Schwarz mechanism. The general form for the low-energy 
charge may be written 

X = mTn + nA + pL, (2.8) 

where Tr is right-handed isospin; A is an axionic PQ symmetry; and L is lepton number, with 
the charge assignments given in the Table. The coefficients m, n and p are subject to (but 
not over-constrained by) several anomaly cancellation conditions [16]. In particular B and 
L violating interactions can be forbidden up to and including dimension six for some choices 
of these symmetries in the supersymmetric limit, as can the /x-terms of the superpotential - 
W ~ UlLH and W ~ fiHH - if n / 0. 
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Very light and weakly coupled gauge bosons from extra-dimensional models 

For the phenomenological discussions of later sections we consider gauge bosons in the MeV 
to TeV mass range, whose direct couplings to Standard Model fermions are much smaller 
than those arising within the Standard Model itself. This section and the next one describe 
several way that very light and weakly coupled bosons can arise from reasonable UV physics. 

Extra-dimensional supergravity provides a simple way to obtain very light gauge bosons 
that are very weakly coupled. A concrete example is six-dimensional chiral gauged super- 
gravity [40], for which the bosonic part of the gravity multiplet contains the metric, (7a/jv, a 
Kalb-Ramond 2-form potential, -B^/jv, and a scalar, (p. Because it is chiral this supergravity 
potentially has anomalies, whose cancellation imposes demands on the matter content. In six 
dimensions Green-Schwarz anomaly cancellation is not automatic, because cancellation of the 
pure gravitational anomalies requires the existence of a specific number of gauge multiplets 
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[41]. Given these multiplets, mixed gauge-gravity anomalies can be cancelled through the 
Green-Schwarz mechanism using the couplings of the field B^n- 

The resulting supergravity admits simple solutions for which the extra dimensions are 
a sphere [42], whose moduli can be stabilized by a combination of background fluxes in the 
extra dimensions [15], and branes coupling to the 6D dilaton [43, 21]. An important feature 
of this stabilization is that the value of the dilaton field becomes related by the field equations 
to the size of the extra dimensions: 

(2.9) 



(M6r)2 



where Mq denotes the 6D Planck scale. This ensures these models are a rich source of U{1) 
gauge bosons, some of whom can have massless modes that survive to low energies below 
the Kaluza-Klein scale. Some of these gauge modes also naturally acquire masses through 
the Stiickelberg mechanism [15] (with the Stiickelberg field arising as a component of the 
Kalb-Ramond field, -Ba/jv)- 

Besides having light gauge bosons, these models also naturally furnish them with very 
small coupling constants. This is because the loop-counting parameter for all bulk interactions 
turns out to be the value of the 6D dilaton, (j), with ~ e"^. But modulus stabilization, 
eq. (2.9), ensures that this coupling can be extremely small because it scales inversely with 
the size of the extra dimensions (measured in 6D Planck units). 

Very light and weakly coupled gauge bosons from low^-energy string vacua 

A related mechanism also often arises in low-scale string models. In early heterotic models the 
role of the Goldstone boson is played by a member of the dilaton super- multiplet: a ~ ImS* 
[35] , while in later Type I and Type II models it is twisted closed string multiplets that instead 
play this role [44, 16]. Although the universal couplings of the dilaton restrict the kinds of 
symmetries that can arise in heterotic constructions of this type, the same is not true for 
Type I and II models. 

There is a simple reason why additional U{1) gauge groups often arise. The basic building 
blocks for constructing models of particle physics in type IIB and IIA string theory are D- 
branes. Generically, the gauge group associated with a stack of N D-branes is U{N), but the 
Standard Model gauge group involves special unitary groups, SU (3) x SU (2) x [/(I). Typical 
GUT models also involve special unitary groups, like 5C/(5), SU{3) x SU{2) x SU{2) x U{1) 
(Left-Right symmetric models) or SU{4:) x SU{2) x SU{2) (Pati-Salam models). It is the 
additional U{l)s that distinguish the Standard Model SU{N) factors from the U{N) factors 
arising from the D-branes, that give new low-energy gauge symmetries. 

Furthermore, anomaly cancellation in string theory typically demands the presence of 
additional D-brane stacks, in addition to those providing the Standard Model gauge group 
factors. These stacks also lead to extra U{l)s under which Standard Model particles are 
charged. Extra U{l)s also appear naturally in F-theory models (for a recent discussion 
see [45]). In many concrete examples these additional gauge fields correspond to U{1)b or 
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U{1)b-l, hence can be relevant for the stabihty of the proton [46, 16] (see also [47] for a 
recent discussion). 

Masses and couplings 

For string vacua the masses and couplings of any gauged U{l)s can be computed, as we now 
briefly describe. 

Consider first the U{l)s associated with the same stack of D-branes as gives rise to the 
Standard Model gauge group. As discussed earlier, such gauge bosons often acquire masses 
from the Stiickelberg mechanism. The size of the mass generated in this way is the string 
scale when the ^7(1) is anomalous [48, 47], but it is the smaller Kaluza-Klein scale for non- 
anomalous U{l)s. 

For models with the compactification volume not too much larger than the string scale, 
these U{1) gauge bosons are very heavy. On the other hand, for large- volume models the string 
scale can be quite low, leading to additional U{l)s potentially as light as the TeV scale. The 
latter can have interesting low energy phenomenology (see for instance [46, 49, 50, 51, 52]). 
In these models the strength of the gauge coupling for the additional U{l)s is roughly the 
same as for the Standard Model gauge couplings (evaluated at the string scale), because both 
have the same origin: the world-volume theory of the stack. Hence they cannot be extremely 
small. 

The masses and couplings of the extra ^7(1) gauge bosons vary more widely when they 
arise from D-brane stacks whose SU (N) factors are not part of the Standard Model gauge 
group. For instance, the case of additional U{l)s associated with D7 branes wrapping bulk 
four cycles of the compactification is discussed in detail in [19]. The value of the gauge 
coupling in this case is inversely proportional to the volume (in string units) of the cycle, S, 
that the D7 brane wraps, 

i-.-^ (2A0) 

In the context of the large volume scenario (LVS) of modulus stabilization [20], the size of the 
bulk cycle associated with the overall volume of compactification can easily be approximately 
Vs ^ 10^ in string units, set by the requirement that one generate TeV-scale soft terms. Thus 
one can obtain gauge couplings as low as g < 2 x 10~^ [19, 53] (couplings larger than this can 
be obtained if the D7 brane wraps a cycle different from the one associated with the overall 
volume). With couplings this small, the gauge boson mass can be Mx ~ gv < 100 MeV even 
if is a TeV. 



3. Gauge boson properties 

With the above motivation, our goal in the remainder of the paper is to work out various 
constraints on the parameters of a massive (yet comparatively light) gauge boson, the X 
boson, that couples to a new U{l)x symmetry. Since the lowest dimension interactions 
dominate in principle at low energies, we include in our analysis all of the dimensionless 
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couplings that such a boson could have with Standard Model particles: i.e. both direct 
fermion-gauge couplings and gauge kinetic mixing. We see how these are constrained by 
present data as a function of the gauge boson mass. 

More specifically, we consider an effective lagrangian density below the supersymmetry 
breaking scale of the form 

C = CsM + Cx + -^mix (3-1) 

where Csm is the usual Standard Model lagrangian; Cx describes the X boson, including 
its couplings to the SM fermions; and £mix is the kinetic-mixing interaction between the X 
boson and that of the SM gauge factor U{1)y [4]. Explicitly 

1 Tn'^ 

£^ = --X^,X^--^X^X^ + iJ^X^, (3.2) 

where X^i^ := d^X,^ — dyX^j^ is the curl of the appropriate gauge potential, X^, and is the 
current for the U{l)x gauge symmetry involving the SM fermions. Similarly, £mix has the 
form, 

Crai. = XB^uX^''' (3.3) 

where is the SM gauge boson for the gauge factor C/(l)y. 

The analysis we provide complements and extends earlier studies of extra gauge boson 
phenomenology. In the lower part of the mass range we may compare with [5], who some 
time ago considered the special cases X = B — L and x = 0- Contact is also possible in this 
mass range with more recent Dark Matter models [9, 2, 3] in the absence of direct matter 
couplings, QxJ^'' = 0. At masses much lower than those considered here other constraints on 
kinetic mixing have also been studied, from the cosmic microwave background [11], and from 
the absence of new long-range forces [31] or milli-charged particles [4, 7]. 

There is also a broad literature on the phenomenology of gauge bosons at the upper end 
of the mass range, largely done in the context of a Z' field and often motivated by GUTs [38, 
54, 55, 56, 57, 58]. Until recently, most did not include the kinetic mixing term. Constraints 
including kinetic mixing arising from precision electroweak experiments are considered in 
[59, 60, 61]; more recent bounds are found in [62, 63, 64, 8]. Many of these analyses overlap 
parts of our parameter space. For instance Z' searches, such as [65], give bounds on the mass 
of the Z' that apply in the regime that the couplings to fermions are identical to that of the 
Z. Others [6, 66] derive bounds for a Z' coupled only to baryon number. 

One difference between the models examined here and those usually considered for Z' 
phenomenology at the weak scale, such as those of ref. [67], is the absence in C of mixing 
between the X and the Z bosons in the mass matrix [i.e. a term of the form £mix = 
6m'^ Z^X'^). We do not consider this type of mixing because we imagine the models of 
interest here to break the X symmetry with a SM singlet. Notice that because the SM Higgs 
is uncharged under the X symmetry, the strong bounds as found, for example, in [68] don't 
apply. 
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3.1 The mixed lagrangian 

In this section we diagonalize the gauge boson kinetic mixing terms (and SM mass terms) 
and identify the physical combination of parameters relevant for phenomenology within the 
accuracy to which we work. Our goal in so doing is to follow ref. [69, 70] and identify once 
and for all how the gauge boson mixing contributes to fermion couplings and to oblique 
parameters [71] modified by the gauge-boson mixing. This allows an efficient identification 
of how observables depend on the mixing parameters. 

We begin by writing the lagrangian of interest more explicitly, after spontaneous symme- 
try breaking. Because it is the Z and photon that potentially mix with the X boson, we also 
focus on these sectors of the SM lagrangian. In order to distinguish the fields before and after 
mixing, where appropriate we denote the still-mixed fields with carets, e.g. X^, reserving 
variables like for the final, diagonalized fields. 

With this notation, the lagrangian of interest is 



where 



with 



•^gauge — •^kin ~l~ ■^mass ('^•^) 



Ain = —^W',.Wr - -^B.uB^" - -^X.^X^'^ + \b,,X^'' (3.6) 
/:mass = -\ (m3#;^ " mo^^) (iTi^W^ - nioB^') - X^X'' , (3.7) 



and 



/ 

Cint = ^ E {^2 {fl^Ts pj) #3 + \JY {Yfp, + yj«7«) /] (3.9) 
/ 

+5x W i^fLlL + X/fi7«) /] X^] . (3.10) 



Here Tsj, Yy^ and Y/-^ denote the usual SM charge assignments, while Xfi^ and are the 
fermion charges under the new U{l)x symmetry. The SM masses, ma and mo, are defined as 
usual [29] in terms of the standard model gauge couplings, gi and (721 and the Higgs VEV, v: 
m3 = 2 g2V and mo = 5 giv. 7^ and 7^ are the usual left- and right-handed Dirac projectors. 
Defining the gauge-field-valued vector V to be 

Y= B , (3.11) 
X 
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the above lagrangian can be written in matrix form 



-gauge 



int 



A H'^ 2 



where 



K 



1 
1 -X 

0-x 1 



and M 



m\ —m^mQ 
-msmo rriQ 
ml 



and 











= E 









_ 92 [fli^Tafji^f] 
91 [/7m (^/i7i + YfRlu) f] 

9X [/7m i^fLlL + ^/K7ii) /] 



(3.12) 



(3.13) 



(3.14) 



The ofF-diagonal elements of M ensure it has a zero eigenvalue, and the condition that the 
matrix K be positive definite requires < 1. 

3.2 Physical couplings 

In order to put this lagrangian into a more useful form we must diagonalize the kinetic 
and mass terms, and then eliminate the SM electroweak parameters in terms of physically 
measured input quantities like the Z mass, Mz, the fine-structure constant, a = e^/Air, and 
Fermi's constant, Gp, as measured in muon decay. 

The diagonalization is performed explicitly in the Appendix, leading to the diagonalized 
form 



4 i^'' 

where the physical masses are 

^2 



1 



(3.15) 



Mi 



and 



Mi 



^1 + s^sh^ry + r^ch^ry + 'dx^ + s^sh^ry + r^ch^ry)^ — 4r^ch^r/^ (3.16) 
^1 + s^sh^jy + richer] — 'dx\J {)- -\- s^sh^ry + r^ch^ry)^ — 4r^ch^T/^ . (3.17) 



In these expressions m| := \ (^i + 52) 

92 

\/9i + 92 
X 



Cy^ '. — COS '.- 



and Sv 



sm( 



91 



while 



shr] := sinhjy := 



X 



and chjy := coshry := 



V9i + 92 
1 



Finally, the quantities and i?x are defined by 

rx ■= and tfx '■= 

ruz 



+1 if rx>l 
-1 if Vx <l 



(3.18) 
(3.19) 

(3.20) 
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which ensures Mz mz and Mx — ?■ mx as — )• 0. 

The currents in the physical basis are similarly read off as 



'11 ■- 



where 



7^ 



-J^s^ + (-J^Svi.sh?7 + J^Cwsh.rj + J^chrj) 
-JT.s 



(3.21) 



E 



Jill [73/7L - Q/si] / 
ie flfiQff 

and e := g2Sw = giCw, ez := e/{swCw) and Q/ = T3/ + 1/^ = Y/^. Finally, := cos^ and 
■s^ := sin^ with the angle ^ given by 

-2s,^sh?7 



^ 1 — s^sh'^r/ — r^ch'^T/ 
Writing the resulting lagrangian as 

C = CsM + ^C,SM + C,X ) 



(3.22) 



(3.23) 



shows that the X boson has two kinds of physical implications: (i) direct new couplings 
between the X boson and SM particles; {ii) modifications (due to mixing) of the couplings 
among the SM particles themselves. 



Modification of SM couplings 

The modification to the SM self-couplings caused by Z — X mixing are given by 



5Cs 



m' 



/ 



(3.24) 



with [70] z := (M| - m|)/m| and 

59fLiR) = (c^ - +5? {sh.r]Sw{QfCw - 5/L(fl)) + ch?? 



(3.25) 



The last step before comparing these expressions with observations is to eliminate the 
parameters and niz (the second of which enters the interactions through r^) from the 
lagrangian in favour of a physically defined weak mixing angle, Sw, and the physical mass, 
Mz ■ This process reveals the physical combination of new-physics parameters that is relevant 
to observables, and thereby provides a derivation [70] of the X-boson contributions to the 
oblique electroweak parameters [71]. 

To this end define the physical weak mixing angle, Sw, so that the Fermi constant, Gp, 
measured in muon decay is given by the SM formula, 

(3.26) 



72 



8s2 M2 ' 
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But this can be compared with the tree-level calculation of the Fermi constant obtained from 
1^-exchange using the above lagrangian, giving (see Appendix) 



-2 _ 2 



1+ ^''^ 



(3.27) 



to linear order in z (which we assume is small — as is justified shortly by the phenomenological 
bounds). 

Eliminating Sw in favour of Sw in the fermionic weak interactions introduces a further 
shift in these couplings, leading to our final form for the neutral-current lagrangian: 

Cr,c = iezY^lr [(ff/f + A<7/0 iL + + A<7/«) 7«] / , (3.28) 
/ 

where := e/s^Cw and 



The SM couplings are (as usual) gj^' := T^j — QfS% and (7|" := —Qfs\,, while the oblique 
parameters [71] 5, T and C/ are given by 

aS = aU = Q, (3.30) 

and 

aT = -z. (3.31) 

Direct X-boson couplings 

The terms explicitly involving the X boson similarly are 

1 Af ^ — 

= —^X^.X^'' - -f X^X^ + ^ E (^/^7^ + ^/h7«) fX^ , 

f 

with 

kfL{R) = eg chr/5(^Xj^(«) + shr/ £;(Q/4- - ^/^(h)) - e^gf^j^) . (3.32) 



We are now in a position to compute how observables depend on the underlying param- 
eters, and so bound their size. When doing so we follow [70] and work to linear order in the 
deviations, A(7j^(^), of the SM couplings, since we know these are observationally constrained 
to be small. 
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4. High-energy constraints 

This section considers the constraints on the X boson coming from its influence on various 
precision electroweak observables measured at high-energy colliders. There are two main 
types of observables to consider: those that test the changes that X-boson mixing induces 
in SM couplings; and those sensitive to the direct couplings of the X boson to SM fermions. 
We consider each type in turn. 

We begin with two well-measured observables that are sensitive only to changes to the 
SM self-couplings: the W boson mass, Mw, and the Z-boson branching fraction into leptons, 
T{Z — )• Later subsections then consider reactions to which direct X exchange can 

contribute, such as the cross section, (Jres(e^e~ — >• /i), for electron-positron annihilation into 
hadrons evaluated at the Z resonance. 

Consistency limits on accessible parameter space 

Since the SM is in such good agreement with experiment [72], it is useful to linearize cor- 
rections to the SM parameters as we have done in the previous section. To be consistent, 
we limit ourselves to considering the subset of parameter space which is consistent with this 
linearization procedure. In practice, we require that the following two conditions of z be 
satisfied: 

1. z must be real (see the discussion in the Appendix), which amounts to demanding that 
it is obtained by a physically allowed choice for the initial parameters rux and x- This 
implies that 

A|-/?|4sh%>0, (4.1) 
where iS.x is defined in eq. (A. 28). This simplifies to 

I Ax - k\> v^k(k + 1) (4.2) 

where 

K:=sls\i^'n. (4.3) 

2. z must be small: z <C 1. To quantify this statement, we assume that z (or, equivalently, 
aT) will be at most within 2(t from its global fit value [70]: 

|z|< 0.014. (4.4) 

This bound has been considered in [73] in the context of hidden sector dark matter 
models. 

In figure 3, we show the regions in the M^— sh r/ parameter space that are excluded by each 
of these bounds. From this, we see that the first condition is dominant when sh?] < 3 x 10~^, 
whereas for greater values of shr/, it is the second condition that is dominant. 
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Plot of the Consistency Bounds on z 
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Figure 3: Plot of the bounds on z as a function of Mx and sliry. Tlie blue crosses limit the region in 
which z is real, and the red squares limit the region in which z <^\. The hatched regions are excluded. 

4.1 Effects due to modified W,Z couplings 

We start with several examples of constraints that probe the induced changes to the SM 
self-couplings. 

The W mass 

Mixing with the X boson modifies the SM prediction for the W mass due to its contribution 
to the electroweak oblique parameter T, as follows [70, 71]: 



where {M'^)sm is the full SM prediction, including radiative corrections: {M'^)sm = -^1(1 — 
loops. Because both the SM radiative corrections and the oblique corrections are known 
to be small, we can neglect their product in the above expression. 

At this point one might ask why bother examine the W mass correction separately, since 
the W mass is one of the observables included in the global fits to oblique parameters, and we 
have already assumed that z must be small enough to ensure that the oblique parameter T 




(4.5) 



(4.6) 



(4.7) 
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Figure 4: Plot of the EWWG bound on the S and T obhque parameters, showing how T is more 
tightly constrained given prior knowledge that S = 0. 



lies within its 2-a range obtained from global electroweak fits (as in Figure E.2 of [72]). The 
reason we re-examine the W mass is that it leads to a slightly stronger constraint, because 
the mixing between the Z and X bosons does not contribute to the S parameter, and this 
prior information leads to a slightly stronger limit on T (as is shown in Figure 4). 

Using the result, eq. (3.31), aT = —z together with eq. (A. 27) for z as a function of r] 
and Mx gives the desired expression for AM^^ as a function of r/ and Mx- In the limit when 
the Z and X masses are very close to one another — i.e. when is such that the equality 
in eq. (4.2) holds — the expression for z becomes 



K - A 



z = = —"OxSwlvl + ^wV"^ + ^iv^) (near-degenerate Z and X masses) , (4.^ 



l + Ax 
and so 

|AM;^| ~ Mzc 



2 



2.75 GeV ( ^ ) . (4.9) 



Moving away from degeneracy, we find the expression for z can be simplified as follows: 



1- Wl-^ 



1 + A;, 1 



+ 0(^2) = + 0(r,^) , (4.10) 



.Rj,\ 2Ax ' ' l-Rl 
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Plot of the Constraint on the W Mass 
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Figure 5: Constraint obtained from limiting the influence of kinetic mixing on the SM value of the 
W mass. The hatched regions are excluded. 



where Rx ■= Mx/Mz {c.f. eq. (A. 22)). So when Mx and Mz are very different, 

AM,, . ,(1^ (4^) . .10 X (^) GeV' . (4.n, 

The large- limit of eq. (4.11) agrees with the result given in [64], which finds 




(4.12) 



The experimental agreement of the measured W mass with the SM prediction implies 
AM^r < 0.05 GeV [74] (2cr uncertainty), and the constraint this imposes on shr/ as a function 
of Mx is shown in Figure 5. Several points about the comparison given in the figure are of 
note: 

• The W^-mass bound on r] is model-independent inasmuch as it relies only on the kinetic 
mixing and does not depend at all on the fermion quantum numbers to which X couples; 

• The strongest constraints on r] occur for Mx nearest to the Z pole, where |f/polc| ^ 
1.8 X 10-3; 

• When Mx ^ Mz, the bound on rj becomes approximately Mj^-independent: \7]\ < 
6.2 X lO"^. This behaviour is also visible in the analytic expression, eq. (4.11); 
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Plot of the Z Decay Bound: shi) = 0.001 



Plot of the Z Decay Bound: sh^ = 0.01 
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Plot of the Z Decay Bound: shi] = 0.1 
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Plot of the Z Decay Bound: shri = 1 
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Figure 6: The constraint arising from Z — ^ £~^£~ decay on the couphng ax = g\/'^T^ as a function of 
Mx, for various values of sh??. The parameters agreeing with the positive bound (AT = +Aroxp) are 
marked with blue crosses, while those agreeing with the negative bound (AT = — AFexp) are marked 
with red squares. The plot assumes a coupling X^^ = X£„ — —1, such as would be true if X = B — L. 
Hatched regions are excluded. 



When Mx ^ Mz, the W mass bounds the ratio ij/Mx- giving Mx/i] 1-5 TeV. 



Z Decay 

The Z decay rate has been measured with great accuracy at LEP and SLC (for details 
regarding their analysis, see [72]). The experimental value [74] for the decay Z — )■ i^£^, 
where i can be any of the charged leptons, is T^+g- = 83.984 =b 0.086 MeV (1 a), and agrees 
well with the SM result [74] 83.988 it 0.016 MeV. The modified Z-fermion couplings change 
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Figure 7: Plot of the constraint arising from considering Z decay into leptons in the hmit where 
Mx ^ ^z- The upper bound (AF = +Arexp) is marked with blue crosses; the lower bound (AF = 
— AFexp) is marked with red squares. Hatched regions are excluded. 



the tree-level decay rate, 
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i+i- 



247r 



(4.13) 



where the couplings Qij 
The deviation from the SM prediction therefore is 



gfj' + A.g£j (with I = L,R) are defined by the interaction (3.28). 
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[2<7,r + Ag,,]A5,,. 



(4.14) 



Notice that this vanishes if ^gii = or when Ag^j = —2gf^'. It can therefore happen that 
Ar£+£- vanishes for two separate regions as one varies through parameter space. 
To obtain bounds on r/ and Mx we use eq. (3.29) to eliminate ls.gfi^i^j^-^^ giving 
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(4.15) 



Here Sj^ + s^, is the SM contribution, gfj^', where 6il denotes a Kronecker delta function. 
Requiring AT£+£- to be smaller than the experimental (2 a) experimental error gives the 
desired bound on the parameters gx, rj and Mx- Figure 6 shows the excluded values in the 
ctx = 5x/47r vs Mx plane, with the leptonic X-boson charge assumed to be = X^j^ = — 1 



(such as would apply if X = B — L). Each panel of the figure corresponds to a different choice 
for sh?7. For the panel in which shry = 1 bounds at lower mass scales than roughly 385 GeV 
are not plotted, since these would conflict with & z = —aT satisfying the global electroweak 
fit, as outlined in figure 3. 

In order to understand the features present in the plots it is useful to consider the small-ry 
limit of z and ^. As discussed above for the W mass bound, the small-77 limit when Mx and 
Mz are very similar or very different must be considered separately. The expressions when 
Mx and Mz are very different are 



As might be expected, all terms in Agij are suppressed by a factor of — M'^/M^ and 

so go to zero when Mx S> Mz- In the opposite limit, Rx — 0, Ag^j ~ Xi^rjSwigx / e-z) + 
s'wC^ / {c^ — •Sw); which can pass through zero (if X^jT] < 0) when \X£j\{gx/ez) — 0{r]). 
Several features of these plots should be highlighted: 

• The best bounds come for Mx — Mz, even for small couplings gx, because in this 
limit the Z — X mixing parameter ^ becomes maximal (tan 2,^ — )• 00), leading to strong 
constraints. 

• For a similar reason, once t] is sufficiently large (shr^ ~ 0.1 — see also figure 7) the 
regime of vanishingly small Ox remains excluded because Agr^^^^) is dominated by the 
oblique corrections to the weak mixing angle. 

• For Mx ^ Mz the excluded area approaches a straight line, corresponding to a bound 
on the ratio gx/M^, as expected from the form of Ag^^^^j^y 

• The graph is more intricate for Mx <C Mz, with slivers of allowed parameter space 
emerging for a narrow, r/-dependent but Mx-independent, value of ax- This happens 
(for sufficiently large r]) because AF = is a multiple-valued condition on the parame- 
ters, as discussed above. 

Figure 7 provides a view of the bounds taken on a different slice through the three- 
dimensional parameter space (77, ax, Mx)- This figure plots the constraints on ax vs shr/, 
in the regime where Mx <^ Mz, showing how a wider range of ax is allowed as shr/ shrinks. 
Note that bounds are only shown for the region where z <^ 1- 

4.2 Processes involving X-boson exchange 

In this section we consider precision electroweak observables, like the resonant cross section 
for e"*"e~ — )■ hadrons, that receive direct contributions from X-boson exchange, in addition to 
the modifications to SM Z-boson couplings. 
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Figure 8: Relevant tree-level Feynman diagrams corresponding to electron-positron annihilation into 
fermion-antifermion pairs. 



The annihilation cross section 

We again proceed by computing the leading change to the tree-level cross section for e~^e~ — >• 
// at leading order in the new interactions. Interference terms between SM loops and X- 
boson contributions may be neglected under the assumption that their product is negligible 
[70]. The relevant Feynman diagrams are shown in Figure 8, where the exchanged boson is 
either a photon, Z or X boson. 

Neglecting fermion masses the relevant spin-averaged squared matrix element for this 
process is (see, e.g. [29] for a treatment of SM scatterings using similar conventions) 

\Y,\^\' = Nc \(\A,, {s)\^ + \A^^ {s)\') + (\A,^ (s)l^ + \A^, {s)\') , (4.17) 



4 

where s, t and u are the usual Mandelstam variables and 

'^"^'>'-= s + S-MI + iV^M, ^ S-MI + iV^M^ • ^ 
The total unpolarized cross section that follows from this is 

a (e+e- ^ //) = ^ (\A,^\^ + \A^^\' + \A,^\' + \A^,\^) . (4.19) 

The couplings gfi and kf, in these expressions are defined in terms of ry, gx and Mx by 
eqs. (3.29) and (3.32). The quantities F^ and Tx are only important near resonance, and 
denote the full decay widths for the Z and X boson, respectively: 

r. = 4^ E [9%+9UN. (4.20) 

2mf<Mz 

^- = ^ E [k). + k),]N,, (4.21) 

2mf<Mx 



where Nc is the colour degeneracy for fermion /. 
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Plot of ai„j(s = M;) Bound: shv] = 0.001 
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Figure 9: The constraint obtained from auad evaluated for ^/s = Mz, as a bound in the ax — 
plane for various values of shr;. Blue crosses (red squares) indicate parameters where predictions differ 
by 2 cr from experiment on the upper (lower) side. The hatched regions are excluded, while diagonal 
shading indicates a region excluded by global fits to oblique parameters. 
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The Hadronic Cross Section at the Z Pole 

Summing the above over all quarks lighter than Mz and evaluating at \fs = Mz gives the 
leading correction to the resonant cross section into hadrons, ahad (s = which is well- 

measured to be 41.541 it 0.037 nb [74]. Requiring the deviation from the SM to be smaller 
than the 2a error gives the desired constraints. Figure 9 shows a number of exclusion limits 
for the coupling ax vs the X-boson mass (for Xj^ = Xfp, = {B — L)f, and Mx in the range 
of 10 — 10^ GeV), with each panel corresponding to a different choice for rj. 

These plots reflect several features seen in the analytic expressions for the couplings: 

• For Mx <C Mz and when rj is small enough, the mass dependence of the bound on ax 
completely drops out, leaving ax 10^^ in this limit. For larger rj small values of ax 
can still be ruled out because the contributions of mixing are already too large. This 
mixing also ensures that the region near Mx = Mz tends to give the strongest bounds. 

• The regime Mx ^ Mz similarly constrains only the combination M'^/ax ^ 800 GeV 
(when rj is small). 

• For T] not too small and Mx smaller than Mz , figure 9 shows a window of unconstrained 
couplings, for the same kinds of reasons discussed above for T£+£- . 

Figure 10 shows a sample slice of the constraint region in the ax vs shr^ plane, in the 
limit Mx ^ Mz- Once again, bounds are not plotted within regions of parameter space for 
which z is not <^ 1. This plot shows that the smallest r/ for which small ax can be ruled out 
is shi] ^ 0.06. Once i] is larger than this, mixing rules out the X boson even with arbitrarily 
small gauge couplings. 

5. Constraints at intermediate energies 

Better constraints on lower-mass X bosons can be obtained from low-energy scattering of 
muon neutrinos with electrons and nuclei. The purpose of this section is to quantify these 
bounds by identifying how the cross section depends on the parameters gx, rj and Mx- We 
consider electron and nuclear scattering in turn. 

5.1 Neutrino-electron scattering 

The Feynman graphs relevant for I'^e" scattering are those of Fig. 8, with three changes: 
(i) the gauge bosons are exchanged in the t-channel rather than s-channel; (ii) there is no 
photon-exchange graph and (in) omission of right-handed neutrino polarizations. 
Crossing to t-channel can be obtained by performing the following substitution 




(5.1) 
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Figure 10: Plot of the constraint from aiiad {s — Mf) in tlie region where Mz- The parameters 

agreeing with the positive bound are marked with blue crosses, while those agreeing with the negative 
bound are marked with red squares. The hatched regions are excluded. 



among the Mandelstam variables in the invariant amplitude ^ X] l-^l^- With these replace- 
ments, the differential cross section for the process f^e~ — >• t'^e" is 



da 

dt 



87rs2 



A,,{t)\^s^ + \A^,{tr{s + ty 



where 



■t-M2 t-Ml' 



(5.2) 



(5.3) 



In the rest frame of the initial electron s ~ 2meEi, and t ~ —2ymeE^, where Ei, is the 
incoming neutrino energy and y is the fractional neutrino energy loss, y := eI /Ey where eI 
is the energy of the outgoing electron. (In such experiments [75, 76], E^, eI ~ 1 — 10 GeV 
so ratios of the form me/Ey and rrie/El can be neglected.) In terms of these new variables 
the differential cross section is 



da 

Ay 



nieEy 



A^MyEu)] + AnMyEy)] {l-yf 



(5.4) 



The cross section for anti-neutrino scattering is easily found from the above by inter- 
changing All ^ 
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Special case: Low-energy limit with rj = 

One case of practical interest is when the boson masses, Mz and Mx, are much greater than 
the invariant energy exchange in the process of interest {i.e. ^y\t\ <C Mx,Mz). When this 
holds the amplitudes, Aj^, can be simplified to 



M| Ml 



(5.5) 



allowing the effects of X-boson exchange be interpreted as an effective shift in the electron's 
electroweak couplings. For E'jy ~ 1 GeV and y order unity this approximation remains good 
down to Mx ^ 30 MeV. 

The resulting cross section is particularly simple in the case of no kinetic mixing, for 
which we can substitute the SM values gei = + and g^i^ = ^ and the X-boson 

couplings fee/ = 9xXei and k^j = Qx^^j, and obtain 

usmg the SM result 2\/2G^ ~ e|/2M|. We see that the X-boson contribution can be 
regarded as an additional contribution to in this limit. This is convenient because it 
allows the simple use of constraints on s\r to directly constrain the ratio g\/Mx- 

The bounds are usually taken from the following ratio [76] of total cross sections, 

R := ''t'^l ^ -'"2 ■ (5-7) 



Given the differential cross section 

[9L+9ln{l-y?] , (5.8) 



dcr , _ _s 2G|me£'^ r 2 2/1 \2i 



dy ^ ^ ^ ' vr 
the total cross section becomes 

a[v^e ^v^e ) = ( 9eL + ) , (5-9) 

and so 

a {V,e- V,e-) = + gl, ) . (5.10) 



Specializing to SM couplings the result depends only on Sy\/ '. 

2>gl^ + gl^ ^ 3-124 + 164 ^ l + K + ^2 

9l, + ^9lu 1-4.2^ + 16.1 1-K + AC2' ^ 



where k = 1 — As^y <ti 1. 
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Plot of the Constraints Associated with Neutrino-Electron Scattering 
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(GeV) 

Figure 11: Bound obtained on ax = Sx/^ti" by limiting the influence of the X boson on the v — 
cross section ratio i?, obtained as a function of for various values of sh -q. The vertical line indicates 
the region ruled out by electroweak oblique fits when 77 = 1. 



Using the experimentallimit [75] As^ = 0.0166 (2a error) with Gj. = 1.1664x10"^ GeV"^ 
[74] to constrain As^ = g\l2\piGpM\ (assuming the choice Xe/Xj^^ = 1, as would be true 
for X = B — L for example) , gives [77] 



9x 



> 4 TeV . 



(5.12) 



General Case: 77 / 

More generally, the couplings kf; also acquire contributions from Z — X mixing even when 
Qx = 0, as the above calculations show. In this case the more general bounds on gx, ij 
and Mx can be extracted by demanding that these contribute within the experimental limit 
AR. Since the experimental limit is often quoted in terms of [75], we translate using 
AR = \dRsM /ds'^ \ As^. In obtaining R, we integrate over y using ^y\t\ ^ M|, but without 
assuming that ^ ^'^x- When evaluating R, we set E^, to a nominal value of 1 GeV. 
Figure 11 shows the resulting bound in the ax — Mx plane, for several choices for r] assuming 
XelX,yL = 1- 

The resulting curves inspire a few comments: 
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For large Mx the bound is independent of r] due to the Mx/Mz suppression of the 
mixing in and A/c. This allows the direct g\lM\ term to dominate. The bounds in 
this regime are relatively strong, and compete with those found in direct searches (e.g., 
by CDF [78] in the case of a SM-like Z'). 

For smaller Mx^ it is the terms in the couplings that are linear in r\ that influence the 
deviation from the 77 = result. To see this, note that 

l^g,^ = -T]—gxX,^ + 0{T]^), (5.13) 

and so there will be a term in |j4j^P that is linear in r] with the parametric dependence 
Qx/Mx- When gx ^ 0, it is this term that is dominant compared to the g\/Mx term 
from X-boson exchange. However, when gx ~ 77, this new term is no longer dominant 
and the bound regresses back to its original slope from the 77 = case at high masses. 



• Once 1 < Mx < 10 MeV and Mx ^ v 1^1; the bound loses its dependence on Mx and 
levels out to some fixed value. This is expected from the form of eq. (5.3). 

• When sh?7 = 1, much of the parameter space is excluded due to the requirement that 
2; <C 1. Therefore, only a small region with > 385 GeV is bounded by electron- 
neutrino scattering in this case. 

5.2 Neutrino-nucleon scattering 

For the bounds from neutrino-nucleon scattering it is worth first recalling how the standard 
analysis is performed. In terms of the neutral-current quark couplings, the quark-level cross 
sections for neutral-current muon-neutrino scattering are 

a (z/^n ^ u^u) = ao (^gl^ + , a {u^d ^ v^d) = ctq (^^L + ^) (5.14) 

a [V^u V^u) = (Jo + gl^ , a {v^d V^d) = ctq + gjj^ 

while those for charged currents are 

a [y^d ijTu) = do and a {v^u — >• /x'^d) = , (5.15) 

where do := 2NcG'],mf,Ey/'K and Nc = 3. 

These show that the quark neutral-current and charged-current cross sections are all 
proportional to one another. The resulting cross section for neutrino-nucleon scattering in 
the deep-inelastic limit is obtained by summing incoherently over the quark contributions, 
giving 

a {v^,N v^X) = ela [v^N ^ X) + e\a [v^N ^u+X) 

a {V^N ^ V^X) = ela {u^N ^ X) + ela {u^N fi' X) , (5.16) 
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Plot of the n-Nucleon Bound: shi^ = 0.1 
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Plot of the n-Nucleon Bound: shii = 1 
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Figure 12: Plot of the constraint from R~ (neutrino-nucleon scattering) assuming X = B — L. Here, 
we plot the bound on ax as a function of Mx for various values of -q. Blue squares (red crosses) indicate 
parameters whose predictions lie 2 a above (below) the central experimental value. The vertical line 
indicates the region excluded by precision oblique fits. 



where 



^l{r) '■— 9ul{r) + 9dL{R) ■ 

The experimental bounds come from the following ratios: 



a {V^N ^ /i+X) 



2 , 2 



(5.17) 



(5.18) 



where r := a {f^N — )■ fJi-^X) ja {v^N — J- ^ X). Most useful is the Paschos-Wolfenstein ratio 



[79], from which the comparatively uncertain ratio r cancels: 



R- 



1 - r a {v^N ^ ^Ji-X) 
Experiments measure the following values [80] 

el = 0.30005 ± 0.00137 
el = 0.03076 ± 0.00110 



(5.19) 



(5.20) 



To constrain the X-boson coupling parameters we work in the regime with y/—t <C M^, 
for which the effects of X-boson mixing and exchange can both be rolled into a set of effective 
neutral-current couplings. The cross sections for quark- level scattering are then given by 
integrating eqs. (5.4) using (5.5), leading to expressions identical with eqs. (5.14) but with 



9ql 



iJql 



9ql9uL + 



Ml 
Ml 



(5.21) 



where q = u,d and I = L,R. Using these in eq. (5.19) gives constraints on e^ = [g^f + 

Figure 12 plots the constraint found by requiring Ae^ < 0.00137, assuming that X = 
B — L. The plots are cut off at low mass where the condition |f/M^| < 0.01 breaks down. 
Notice that for r/ = the bound is similar to that found for neutrino-electron scattering, 
with stronger bounds on ax at smaller Mx- For nontrivial r] the strength ol X — Z mixing 
eventually provides the strongest constraint, leading to strong bounds even for small gx at 
sufficiently low Mx- 

6. Low-energy constraints 

We finally turn to constraints coming from lower-energy processes. 
6.1 Anomalous magnetic moments 

The accuracy of anomalous magnetic moment (AMM) measurements [74] produce a strong 
constraint on the parameters of an extra gauge boson. We consider the bound arising from 
both the electron and muon AMM on the X gauge coupling as a function of the mass Mx, 
for various values of the kinetic mixing parameter sh?]. 

The correction to the AMM of a lepton, £, is given by [3] 



6ap 



An'^Ml 



'^^ klzjl - zf - kj^ [z{l -z){3 + z) + 2(1 - zfmj/M^] 
^ z + (1 - z)2m2/M2 



where the vector and axial couplings to the X boson are of the form 



kiL - ken 



chijgxXiv 
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(6.1) 



(6.2) 
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Plot of the AMM Bounds: shrj = 




Plot of the AMM Bounds: shn = 0.001 
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Figure 13: Plots of the constraint on the gauge coupling ax arising from the electron and muon AMM 
as a function of M^, for various values of shry. The electron AMM bound is marked with blue crosses; 
the muon AMM bound is marked with red squares. The plot assumes a coupling = Xgi^ = — 1, 
such as would be true \i X = B — L. Hatched regions are excluded. 
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There is, however, some subtlety in comparing this shift with experiment [10]: since the 
electron AMM, 6ae, is used to determine the fine-structure constant, a. The best bound on 
X boson couplings therefore comes from the next most precise experiment that measures a, 
and not the errors from the {g — 2) experiments themselves. Following [10] this leads to the 
constraints 6ae < 1.59 x 10"^'' and (5a^ < 7.4 x 10~^, which when compared with the above 
expression gives the bounds shown in Figure 13. These plots reproduce the results found in 
[3] when shr/ = 0. In particular, the Mx values below which any gauge coupling is excluded 
are consistent with the bounds shown in [10]. 

Since these bounds are often considered (e.g. in [10], [9, 2]) in the context of a constraint 
on kinetic mixing, we also plot the constraint on shr/ as a function of the X boson mass, for 
various values of the gauge coupling. This is shown in Figure 14. 

6.2 Upsilon decay 

The bound we present here is an extension of the result found in [2] . By looking at the decay 
rate of the T(3s) bb bound state, researchers from the BABAR collaboration were able to 
place a bound on the occurrence of a particular channel involving a light pseudoscalar Aq 
[81]: 

e+ + e-^T(3s)^7 + Ao^7 + /^^ + Ai"- (6.3) 
Their upper limit on the number of events 

AT = a{e+ + e- ^ T(3s)) x £ x Br(T(3s) ^ 7 + Aq) x Br(Ao ^ /x+ + /x") , (6.4) 

places a bound on the quantity Q := Br(T(3s) — )■ 7 + Aq) x Br(^o — 5- + fJ-~)- 
However, the reaction of interest to us is 

e+ + e- ^7 + X^7 + ;U+ + /i", (6.5) 

which would have an identical signature. So the measured bound can also be reinterpreted 
as applying to the quantity 

- :'(::::-:;;.))' <-) 

The experimental limit [81] Qx < 3 x 10~^ gives the plots found in Figure 15 over the range 
2?Ti^ < Mx < Ecm{= 10.355 GeV). This bound is quite strong, as it eliminates the entire 
region for shry > 0.002 (as is shown in [2]). For smaller shr/, the bound is roughly constant 
when Mx <C Ecm- 

As in the case of the AMM bounds, we also plot the constraint on sh?/ as a function of 
the X boson mass for various values of the gauge coupling — as shown in Figure 16. 



- 39 - 




Figure 14: Plots of the constraint on the kinetic mixing, shry, arising from the electron and muon 
AMM as a function of A/x, for various values of the gauge coupling ax- The electron AMM bound 
is marked with blue crosses; the muon AMM bound is marked with red squares. The plot assumes a 
coupling — Xgjf = — 1, such as would be true ii X = B ~ L. Hatched regions are excluded. 

6.3 Beam-dump experiments 

In the MeV— GeV mass range, small gx and rj are constrained by several beam dump exper- 
iments. These bounds are considered in detail in [2]; we apply a simplified version of their 
analysis here. 

In these experiments, a large number of electrons with initial energy E are collided 
with a fixed target made of either aluminum or tungsten. Many of the resulting collision 
products are absorbed either by the target or by some secondary shielding. (Here, we use t 
to denote the total thickness of both the target and the shielding.) The remaining products 
continue along an evacuated tube to the detector, located at some distance D away from the 
target. For a summary of values for these parameters, see Table 1. 
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Figure 15: Plots of the constraint on the gauge coupUng arising from T(3s) decay as a function 

of AIx, for shry = 0,0.001. The plot assumes a coupling = Xi^ — —1, such as would be true if 
X = B — L. Hatched regions are excluded. 




(GeV) (GeV) 

Figure 16: Plots of the constraint on the kinetic mixing, sh77, arising from T(3s) decay as a function 
of Mx, for = 0, 1 X 10^^. The plot assumes a coupling X^l = Xen = —1, such as would be true if 
X = B — L. Hatched regions are excluded. 

The bound arises from the non-observation of X decay products. The incoming electron 
emits an X boson as bremsstrahlung during photon exchange with the nucleon (N): +N ^ 
e~ + N + X. The X can then decay into either an e^e~ or /-i^/i^ pair. However, a decay 
that occurs too soon is absorbed by the shield while a decay that occurs too late occurs past 
the detector. Therefore, the number of lepton anti-lepton pairs observed at the detector can 
be computed by multiplying the number of X bosons produced, Nx, by the probability for 
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Experiment 


Target 


Ne 


Beam Energy 


t 


D 


E774 


W 


0.52 X 10^° 


275 GeV 


30 cm 


7.25 m 


E141 


w 


2 X lO^^ 


9 GeV 


12 cm 


35 m 


E137 


Al 


1.87 X 10^0 


20 GeV 


200 m 


400 m 



Table 1: Parameter values for the E774, E141, and E137 beam dump experiments. 




Figure 17: The constraint arising from beam dump experiments on the coupling ax — dx/^^ ^ 
function of Mx, for sh?7 = 0, 0.001. The E774 bound is marked with red squares; the E141 bound is 
marked with blue crosses; the E137 bound is marked with black circles. The plot assumes a coupling 
= Xiji = — 1, such as would be true if X = B — L. Hatched regions are excluded. 



the X to decay between z = t to z = D: 

/ 1 

■ e~'/"^ ] . (6.7) 



Here, we write the lab frame decay length as Iq := 'Jct, where 7 = (1 — the relativistic 

time-dilation factor and r is the inverse of the X rest-frame decay rate: r := 1/Tx- 
In estimating the number of X^s produced, we use the following result from [2]: 

Ml ' 



^x-^e/i'^, (6.8) 



,2 o K i\/r„AA2 



where e = x^w and fi ~ 2.5 MeV is an overall factor that contains information regarding 
the details of the nuclear interaction, and is shown in [2] to be roughly constant for Mx 
between 1 and 100 MeV. There is, however, an obstacle in applying this result directly to 
our analysis: it was derived without including any coupling to J^. In order to introduce 
the A;eL(H)"dependence in this expression, we note from [2] that the e-dependence above arises 
from the cross section a{e~j — )• e"X) under the assumption that the electron is massless. 



This means that the left- and right-handed hehcity X — e interactions contribute equaUy to 
the cross section, ahowing the substitution 



with the normahzation chosen so that the above expression reduces to x^c^, in the case where 
XeL{R) = 0, sh?7 <C 1 and Mx <^Mz. 

All in all, we find that the number of X's we expect to observe is given by 

Applying the experimental exclusions [2] Nobs < 17 events (E774), Nobs < 1000 events (E141), 
and Nobs < 10 events (E137) gives the bounds shown in Figure 17. 

The plot for shry = gives good agreement with a similar plot in [3] (in the region over 
which these results overlap). The lower bounds for each experiment are approximately flat 
because, in the region where t <^ D <^ Iq, the fraction of X's that decay is just D/io, which 
gives 



The leading Mj^-dependence then cancels since ~ 1/M^. The upper bound results from 
the situation where the X bosons decay too quickly, and the decay products do not escape 
the shielding. 

We have only included plots for the cases where shr/ = and 0.001 because the bounds 
become too weak to constrain any region of this parameter space whenever sh?] > 0.007. 

An interesting feature of these bounds is that, at any given value of sh?7, the gauge 
coupling can be increased such that the bounds are evaded. This occurs because a stronger 
gauge coupling causes the X bosons to decay within the shielding. Therefore, any bound 
on kinetic mixing which results from these experiments can weaken if the direct coupling of 
electrons to the X is taken to be non-zero. To demonstrate this, consider the bounds shown 
in Figure 18, which plots the bound on kinetic mixing as a function of the X-boson mass, for 
various values of ax- Note that, for ax ^ 1 x 10"^, these bounds are satisfied for all values 
of sh?7 in the relevant mass range. 

6.4 Neutron-nucleus scattering 

Low-energy neutron-nucleus scattering is important because most of the other low-energy 
bounds evaporate if the new boson doesn't couple to leptons (such as if X = B). For 
neutron-nucleus scattering a bound is obtained by considering the effects of the new Yukawa- 
type potential that would arise from a non-zero vector coupling of the X to neutrons. For 
light X bosons this can be seen over the strong nuclear force because it has a longer range, 
and can affect the angular dependence of the differential cross section for elastic scattering. 
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Figure 18: The constraint arising from beam dump experiments on the kinetic mixing shrj as a 
function of Mx, for = 0, 1 x 10~^^, 1 x 10~^^, and 1 x 10~^. The E774 bound is marked with red 
squares; the E141 bound is marked with blue crosses; the E137 bound is marked with black circles. 
The plot assumes a coupling X^,^ = X^^ = ^ 1, such as would be true if X = B — L. Hatched regions 
are excluded. 



da{nN — )• nN)/dO,. This bound is discussed in the context of a scalar boson in [82] and more 
generally in [83]. 

Following these authors we parameterize the differential cross section as 

^ = ^(l + .Eco.«), (6.12) 
where ctq and w are to be taken from experiments. Then an interaction of the form 

^W0= |J — (6.13) 
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leads to a correction to the expected value of w, which is measured experimentally in the 
energy range E ~ 1~10 keV for neutrons scattering with ^^^Pb. 
Agreement with observations leads to the bound [82, 83] 



k 



< 3.4x10"" 



(6.14) 



where knv = kuv + 2fcrfy with kjy ■= \{kfL + kf^)-, as above. Figure 19 shows a plot of this 
bound for the nominal case shr/ = 0. 

For this combination of couplings, an interesting cancellation occurs. For small kinetic 
mixing, the correction AA;j^(fl) has the form 



= r]ecwQf + 0{rf ) , 



(6.15) 



and so the leading correction in rj vanishes for any electrically neutral particle, like a neutron. 
This makes this bound relatively insensitive to changes in kinetic mixing, not varying appre- 
ciably over the range < sh?y < 1. A similar cancellation occurs in the case of nucleosynthesis, 
considered in S6.6. 



Plot of the Neutron-Nucleus Scattering Bound 




10 



J I I I I I I 



10" 



10" 



10" 



10 



-1 



10" 



MJGeV) 



Figure 19: Plot of the constraint on the gauge coupling ax due to neutron-nucleus scattering as a 
function of the X-boson mass Mx- The hatched regions are excluded. 
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6.5 Atomic parity violation 

The Standard Model predicts a low-energy effective coupling between the electron axial cur- 
rent and the vector currents within a given nucleus. The so-called weak charge of a nucleus 
with Z protons and neutrons is defined (up to an overall constant) as the coherent sum of 
the Z-boson vector couplings over the constituents of that nucleus [84]: 

N) := 4 \Z {2guv + Qdv) + N {g^v + 2gdv)\ ■ (6.16) 

where 

._ 9fL+9fR , ._ SfL- 9fR . „s 

9fvf ■- ^ and qja ■- ^ • (b-i/j 

In terms of these the leading parity- violating effective electron-nuclear interaction gener- 
ated by Z boson exchange is 

Ceff = -V2GpgeAQw (eT^Tse) {'^-f"^) . (6.18) 

where ^ is the field describing the nucleus. X-boson exchange adds an additional term to 
this effective lagrangian of the form 

C^ff = (e7.75e) (^7^^^) (6-19) 

where 

Qx := Z {2kuv + kdv) + N {kuv + 2kdv) ■ (6.20) 
Therefore, the total shift in the Qw due to the X boson is 

2y/2 keAQx 



w 



where 



Ml 



(6.21) 



Qf,«(Z, AT) = 4 [Z (2gl^ + ff^) + iV {g^l + 2^,^)] 

= Z{\-\s\)-N . (6.22) 

Notice that the bracketed term in AQm^ goes to as r/ —t- 0, whereas the second term does 
not as long as fexe does not vanish in the same limit. The total effective lagrangian for this 
system can then be written as 

+ Z:^// = ^ (Qr + AQ,.) (e7^75e) (^7^^^^ ) • (6-23) 

It is expected that the second term in eq. (6.21) will be dominant, so it is useful to 
consider the form of kAe/M^ in the limit where Mx ^ Mz and ?] ^ 1: 

keA _ 9xXeA (l + ^ C^^^) (aoA \ 

G,Ml- G,Ml ^ e,''- ^^-'^^ 

Therefore, if = X^g, then the constraint becomes significantly less stringent at low masses 
since, instead of bounding the ratio g^/Ml, it is now the combination gxTj that is bounded. In 
order to emphasize the strength of this bound when 7^ 0, we use the charge assignments 
as shown in Table 2. 
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Plot of the Constraint from the Weak Charge of Cesium 
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Figure 20: Plot of the constraint on the gauge couphng ax due to the weak charge of cesium as a 
function of the X-boson mass Mx, for various values of sh?/. The hatched regions are excluded. 
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If the X boson is light enough the above effective in- 
teraction eventually becomes inaccurate in describing the 
electron- nucleus interactions. In this case, rather than pur- 
suing a detailed analysis of the microscopic lagrangian, we 
follow ref. [84] and introduce a corrective factor K[Mx) to 
account for the non-locality caused by the small mass of the 
X boson. This modifies our expression for ^Qw as follows: 



Table 2: Charge assignments for 
the "right-handed" C/(l). 



9eA 



-1/4) 



Qw Q 



SM 
W 



G f 



Ml 



K{Mx). (6.25) 



In [84] a table is given for K for various values of Mx in the range 0.1 MeV < Mx < 100 
MeV. In order to render the graphs shown here, we have interpolated values of K by doing a 
least squares fit to the values in [84]. 

As with the neutrino-electron scattering bounds, the slope of the bound changes for ?? 7^ 
due to the production of a new dominant term through cancellation with the modified Z- 
fermion coupling. Once again, we exclude the region below 385 GeV for the sh?/ = 1 plot in 
order to avoid conflict with the electroweak oblique fits that require z <C 1. 

Since this bound relies crucially on there being an axial vector coupling to the electron, 
we did not include it when compiling the summary of bounds given in figures in §1. 
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6.6 Primordial nucleosynthesis 

We close with the study of constraints coming from cosmology, which for the mass range of 
interest in this paper consists dominantly of Big Bang Nucleosynthesis. 

Any X bosons light enough to be present in the primordial soup at temperatures below 
T ~ 1 MeV can destroy the success of Big Bang Nucleosynthesis (BBN) if they make up 
a sufficiently large fraction 10%) of the universal energy density, leading to potentially 
strong constraints. In particular, such a boson poses a problem if it is in thermal equilibrium 
at these temperatures. 

Quantitatively, measurements of primordial nuclear abundances forbid the existence of 
the number of additional neutrino species (beyond the usual 3 of the SM) to be [85] SN^, < 1.44 
(at 95% C.L.). But since each boson in equilibrium counts | times more strongly in the 
equilibrium abundance, and since a massive X boson carries 3 independent spin states, the 
corresponding bound on the number, Nx, of new species of spin-1 particles in equilibrium at 
BBN is 

iV^< 0.84. (6.26) 

Even just one additional massive spin-1 boson into relativistic equilibrium is excluded at the 
95% confidence level. 

In a universe containing only the X boson and ordinary SM particles at energies of 
order 1 MeV, this leads to two kinds of constraints: either the X boson's couplings are weak 
enough that it does not ever reach equilibrium; or if the X boson is in equilibrium it must 
be heavy enough (> 1 MeV) to have a Boltzmann-suppressed abundance. Figure 21 sketches 
the regions in the coupling-mass plane that are excluded by these conditions. The vertical 
line corresponds to the situation where abundance is suppressed by Boltzmann factors. 

The constraints on couplings arise only for sufficiently light particles, and express the con- 
dition that the couplings be weak enough to avoid equilibrium, at least up until the freeze-out 
temperature Tp. There are two curves of this type drawn, which differ by whether it is col- 
lision or decay processes that are the dominant equilibration mechanisms. Qualitatively, the 
requirement that reactions like Xv -H- Xv not equilibrate the X bosons leads to a constraint 
on the couplings that is Mx-independent in the limit where Mx <^ Tp, because then the size 
of both the reaction rate and Hubble scale is set by the temperature. The same is not true 
for decay reactions, X ■^Vi', since the rate for this also depends on the X-boson mass. 

A few other comments are appropriate for Figure 21. First, because they are outside the 
main scope of this study, the bounds shown are derived assuming that Mx <C T (rather than 
being evaluated numerically as a function of Mx) and so are drawn only up to the mass range 
within 0.5 MeV of the freeze-out temperature. Second, the resulting expressions depend only 
weakly on rj, showing little difference over the range < shr] < 1. As discussed in earlier 
sections, this is a consequence of the neutrino's electrical neutrality, which ensures that the 
leading small-77 limit of the kinetic mixing first arises at 0{rf') rather than 0{r]). 
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Figure 21: Constraint on the gauge coupling of the X due to its effect on nucleosynthesis, as a 
function of the X-boson mass. The red squares indicate the bound due to Xv^Xv scattering; the 
blue crosses indicate the bound due to X^vV decay. 
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Note Added 

Since posting, we have learned of a beam dump analysis [86] that has enlarged^ the exclusion 
regions discussed in §6.3. 

A. Diagonalizing the gauge action 

This appendix provides the details of the diagonalization of the gauge boson kinetic and mass 
^We thank Johannes Bliimlein for bringing this to our attention 
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mixings. The starting point is eq. (3.12), 

1 



(A.l) 



with K and M given in eqs. (3.13). 



Diagonalization 

We begin by performing the usual weak-mixing rotation to diagonahze the mass term: 



V = i^iV : = 



where 



Cw '■= coson- :- 



92 



1 
and Sw '■= 



smf 



91 



V9i + 9i 



The lagrangian then becomes 

C = 

with new matrices 



k = RlKRi 



1 xsw 
1 -XCw 
XSw -XCw 1 



and M = RlMRi 



mi 

ml 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



where m| := | (^gf + g(|) v^. Under the same transformation the currents become 



J fl + Jfi Cw 



(A.6) 



le 



which defines := e/{swCw) and uses the standard SM relations g2Sw = 9iCw '■- 
Qf = T3f + Yf, = Yf,. 

The kinetic term is diagonalized by letting 



e and 





'l 





-Sw shr] 




' Z 


V := LV := 





1 


Cw shr] 




A 










chrj 




X 



with 



shrj := sinhry := 



X 



and chry := coshr/ : = 



X^ 



(A.7) 



(A. 
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This gives, by construction 



K = L'^KL = 



1 
1 
1 



and 



M = L^ML 












(A.9) 



(A.IO) 



while the currents become 



7^ 



(A.ll) 



(Notice that L and i?i satisfy LRi = RiL, so it is immaterial whether we first diagonalize 
the SM mass or the kinetic terms.) 

Finally, the mass matrix is diagonalized by letting 



(A.12) 














'z' 


V = R2Y := 





1 







A 









. 




X 



where := cos^ and := sin^ with the angle ^ given by 

—2swsh.ri 



tan 2^ = 
where we define for convenience 

The diagonalized lagrangian then is 



1 — s^sh^ry — r'j^ch^r] ' 



mx 
mz 



4 

where the physical masses are 

Mi = — 
^ 2 



1 



Mi 



1 + s^sh^?7 + r|ch^?7 + i?x\/(l + s^sh^r? + r^ch^ry)^ - Ar'^da^r] 
^1 + s^sh^?7 + r\da^rj — dxsj {)- + s^sh^?? + r^ch^r/)^ — 4r^ch^?7^ 



with i?x defined such that ruz and Mx ^ 'Tlx as 77 ^ 0: 

%1 if rx > 1 



^x := 



-1 if rx < 1 



(A.13) 
(A.14) 

(A.15) 

(A.16) 
(A.17) 

(A.18) 
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The currents in the physical basis are similarly read off as 





' '^<^ii + ( 


—J^Swsh.!] + J^Cwshry + J^cYitj) 






3,= 














[-J^Swsh.r] + J^Civshr/ + J^chrj) 







(A.19) 



Since we are eventually interested in obtaining bounds in terms of the physical masses 
Mz and Mx, it is useful to invert these mass equations to find the input parameters m| and 
as a function of the physical masses and r]. This gives 



Ml 
2ch^r? 



l + Ri + {)xJ{l + Rir -41 + slsh'ii) Rl 



Ml 



2 (1 + sl^h^iT) 

where Rx is used to denote the ratio of the physical masses: 



(A.20) 



l + Rl- -dx^il + Rlf - 4 (1 + slsh^ Rl (A.21 



R. 



Mx 
Mz 



(A.22) 



Also, the sign 'Ox is now +1 for Rx > 1 and —1 for Rx < 1. 

Given this inversion, the angle ^ can now be written as a function of Rx and rj only: 



2s„/sh7/ 



1 - s^sh. 7] - r1{Rx,i])ch 7] 



where 



t&n2^{Rx,v) 



1 + rslsh'r]) 1 + + / (I + RIY -41 + slsh'r]) Rl 



(A.23) 



rxiRx,v) 



ch^T? (l + Rl- ^xJ{l + Rlr - 4 (1 + s^sh\) Rl 



(A.24) 



Physical couplings 

We are now in a position to read off the physical implications of the X boson. That is, we 
may write 

where the modification to the SM self-couplings are given by 



6Cc 



■- m„ 
2 



/ 



(A.25) 
(A.26) 



with [70] 



z{Rx,v) ■-- 



Ml - ml slsl?r] -/^x + ^x \ ^l - Rlslsb?r] 



i + ^x-^xJ^l-RVsU^^ 



(A.27) 
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where 



A,:=^(i?|-1). (A.28) 



(Note that the r/ — )• Umit of z is easily verified by implementing the identity = "dx sj ^x-) 
Given the form of 2, one might worry that, for some choice of the parameters Mx and 
shr/, z would yield a complex value. However, any such choice does not correspond to a choice 
of real values for the original parameters of the lagrangian, m^, w-z, and x- This happens 
because sufficiently large kinetic mixing tends to preclude the existence of mass eigenvalues, 
Mx and M^, that are too close to one another. This is why this region of parameter space is 
excluded from the plots of §4. 

The fermion couplings are similarly 

(^^/LCfl) = (c^ - 1)9/L(i?) +sg (^shrysw'(Q/ci + chr?|^Xji(«)j . (A.29) 

The terms explicitly involving the X boson are 

l^x = -\ X^.X^'' - ^ X^X^ (A.30) 

/ 

with 

^/L(fl) = ( ch-qgxXf^j^) + shr/ A(Q/c^ - 9fL{R)) ) - st^ezgfL(R) ■ (A.31) 

Notice that in this basis X^ does not couple directly to the electroweak gauge bosons at 
tree-level, but has acquired modified fermion couplings due to the mixing. 

Oblique parameters 

The only remaining step is to eliminate parameters like Sw and in the lagrangian in 
favour of a physically defined weak mixing angle, Sw^ and mass M^. This process reveals the 
physical combination of new-physics parameters that is relevant to observables, and thereby 
provides a derivation [70] of the X-boson contributions to the oblique electroweak parameters 
[71]. 

We have already seen how to do this for the Z mass, for which 

~ Mz (1 - I) ■ (^-32) 

For the weak mixing angle it is convenient to define Syy so that the Fermi constant, Gp, 
measured in muon decay is given by the SM formula, 
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But this can be compared with the tree-level calculation of the Fermi constant obtained in 
our model from W^-exchange, 



9l 



to infer 



-s2 ^2 2 2 /-i I \ 

^W^W ~ ^W^W v-*- "T -^j ' 



which, to linear order in z, implies that 



1 + 



{AM) 
(A.35) 

(A.36) 



Eliminating Sw in favour of Sw in the fermionic weak interactions introduces a further 
shift in these couplings, leading to our final form for the neutral-current lagrangian: 



le 



■'NC 



SwCw 



le 



SwCw 



f 

(-|)E{77' 



f 



■■= iezJ27r [{gfr + Agf,)j, + {gf^ + Agf,)j,] f Z, 
f 



(A.37) 



where := e/swCw and 

^9fL(R) 



-,SM 



z 



(A.38) 



where (as usual) gj^ := T^f — QfS%, and := —QfS^. It is assumed throughout that the 
corrections z, Sgfj^f^j^-j, and 

Akfi^i^ji) ■= kfL{R) - £'x^/L(fl) (A.39) 

are small, so that any expression can be linearized in these variables. In particular, this means 
that one can replace hatted electroweak parameters (i.e. s^r, etc..) with unhatted ones in 
our previous expressions to give: 



z{Rx,v) = 



s^sh^r? - A;, + Al - Rlslsh.'^ri 



(A.40) 



l + A^-^x^Al-Rlslsh^T] 
S9MR) = {c^-l)gf^^^)+s^(^shr^s^{Qfcl-gf^^^)) + chr,^^Xf,^^)^ (A.41) 
Afc/^(fi) = (c^chr? - 1) 5xX/i(«) + c^shr] -^{Qfcl - gj^^^^) - s^e^gf^^^) . (A.42) 
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Alternatively, one can use the relationship between z and r] to determine the contribution 
to the oblique parameters [71] 5 = C/ = and aT = —z, where (as usual) a := e^/An. In 
this case, A(7j^(^) can be written as in [70] 

^9M.) = ^ af^in) + (-fi^) Qf + . (A.43) 
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